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Abstract 

We show that given n normahzed intervals on the unit circle, the numbers of visits 

of d random rotations to these intervals have a joint limiting distribution as lengths of 

C/^ I trajectories tend to infinity. If d then tends to infinity, then the numbers of points in 

^H ' different intervals become asymptotically independent unless an arithmetic obstruction 

r~| , arises. This is a generalization of earlier results of J. Marklof. 

-)— > 

The following question arises from two results of Marklof about gap distribution for rota- 
tions. Fix a point ^ in [0, 1) and let Bn = {(},N^/('^-^)]®('^-^') © R c R'^. What what is the 
limiting behavior of the number of points of the form 
> 

(N : fd-i 

Y^ ruiai (mod 1) : m^ G [1, TVVC'i-i)] n Z, 1 ^ i ^ rf - 1 

i=l 






^ '. for ai G [0, 1) that land in (^ - f , ^ + f ) as iV ^ cx)? In ^Ij, J. Marklof showed that 

O' 

leb \ (a,0 e [0, ly-' X [0, 1) : #{m eB^nZ":^ + J^m.a, + m, e {-§, f )} = A [ ^ p('^)(A) 

c^ ' as A^ — !■ oo and found its decay as A — )■ oo. His main tool was the mixing property of 

a diagonal flow on SL((i, R)/SL((i, Z) that had been proved by Moore [3]. In a later note 
Marklof remarked that for one variable (that is, d = 2), a stronger result is true due to a 
Theorem of Shah [B]. Namely, for fixed ^ G [0, 1) \ Q, 

leb {a G [0, l):#{mG {!,..., iV}:e + ma (mod 1) e {-§,§)} = A] ^ P^'\A) . 

This result uses Ratner's Theorem on measures invariant under unipotent flows [4J. We will 
generalize the theorems mentioned above to joint limiting probability distributions for several 
intervals and study their large d limits. 



1 Notation and results 

We will use the following notation. 

• A^, n, and d denote positive integers with d ^ 2; 

• upper indices (usually j) run from 1 to n and lower indices (usually i) run from 1 to d 
unless stated otherwise; 

• m = {rrii, . . . , rrid) is a vector of d integers; 

• (T = {a^, . . . , (7^) is a positive vector {a^ > for all j); 

• (x = {ai,...,ad-i)e{K/ZY-'- 

• ^ = (eS...,r)e(R/Z)"; 

• T = (r^, . . . , r") is a real vector; 

• Pois a denotes the Poisson distribution with parameter a. 

If n = 1, we will write a instead of a^ and cr and similarly for other variables. Let B^ = 
(0, iV^/'^'^"^)]®*^'^"^) © R C R"^ as before. For a measurable set S define random variables 
Xg^:[0,l)^Zby 

We will usually suppress the upper indices on X^^s- For a vector ^ G T", the set H C T*^ 
is the closure of the orbit of rotation by $, on the torus: S = {k^: k G Z}; it is the smallest 
closed Lie subgroup of T" that contains ^. 

Our results for limiting distributions of X^^s are as follows. 

Theorem 1. Fix any absolutely continuous probability measure on [0,1). With notation as 
above, the distribution of 

X^^T,(T = (-^Cl,(rl,Tl+(Tl); • • • ;-^^",(t",t"+o-")) 

has a weak limit as N ^ oo; we denote it by P ^ ,„. 

In other words, the numbers of points in shrinking segments (^■' + ^,^"' + ^-^^), 1 ^ j ^ n, 
with fixed "centers" ^^ have a joint limiting distribution as A^ tends to infinity. The limiting 
distribution depends on d, n, cr, S, and r modulo S. In particular, if S = T", then the 
distribution is independent of r. 

Remark 1. Jens Marklof proved special cases of this theorem. He proved the case n = 1, d = 2 m 
[2] and the case n = 1 and arbitrary d with average over ^ in [1]. 



Theorem 2. Let P ^ ,„ be the distribution from Theorem 1. Then, P ^ ,„ has a weak 
limit as d ^ oo. Furthermore, 

Pi'!l,H,T/H =^ (Pois ai, . . . , Pois an) 

as d ^ oo iff 

{rf T^ + a^) n {t\t^ + a^ ) = whenever ^•' = ^■' ■ 

In effect, this Theorem says that as the number of rotations tends to infinity, the gap 
lengths exhibit random behavior. However, for every finite d and S (with n ^ 2) we have that 

-Pno-Hx/H ^^ dependent. 

Acknowledgements. I would like to thank Ya. Sinai, E. Lindenstrauss, J. Marklof, F. 
Cellarosi, A. Salehi Golsefidy, and Z. Wang for helpful discussions. 

2 Large A^ limit 

Proof of Theorem 1. We reformulate the problem in the language of homogeneous spaces. Let 
L = SL(d, R) K (R'^)®" and let A = SL(rf, Z) k (Z'^)®" C L. Muhiplication law on L is given 
by 

(M, v\..., 'y")(N, w\..., w'') = (MN, v^ + Mw\ ...,«" + Mw""). 

It is well-known that A C L is a non-cocompact lattice. The homogeneous space L/A is a 
bundle over SL(rf, R)/SL(rf, Z) with fiber (T'^)®". 
Given a set of vectors v^, . . . ,v'^ E T'^, let 

Lv = {(1, v\..., v")-i(M, 0, . . . , 0)(1, v^ . . , v") I M e SL{d, R)} C L; 

it is of course isomorphic to SL((i, R). Also define Ly to be the smallest group containing 
Lv that is defined over Q. Dimension of Ly depends on vf If all vectors v^ have rational 
coordinates, then Ly = Ly- Otherwise the fiber over the identity in Ly is the smallest Q- 
vector space containing the identity fiber for Ly. This construction can be carried to other 
points. Finally set Ay = Ly fl A which is a lattice in Ly by construction. 

For our purposes fix v^ = (0, . . . ,0,^-')"^. For this choice of v^ we get the homogeneous 
space Ly/Ay. We have constructed Ly so that n{Ly) = 7i{Ly), where tt: L — > L/A is the 
canonical projection. The structure of this space depends on S = Z^ C T". It is a subbundle 
of L/A: the base is still SL{d, R)/SL((i, Z) but the fiber is S*^ after reordering coordinates. 

We define /^,^: Ly/Ay -> R" by 



where 



meZ''\{0} 

ji xe (o,a) 

Xa{X) = < 

otherwise, 



and 



m = Mm + V. 



It is easily seen that / is Ay invariant and hence well-defined on the quotient. 

We need to show that leb{/^,^ = {A\ . . . , A")} ^ Pi'!l,H,T/H(^^ • • • , ^") as A^ ^ oo. To 
this end we use a theorem of Shah (Theorem 1.4 in [6j). The form we need is the following: 

Theorem 3 (Shah). Let 



Un 



(\ 



\ai 



\ 



ttd-i 1/ 



and $+ 



/e-* 



V 



\ 



^{d~l)t 



Let L he a Lie group, K d L a lattice, ip: SL((i, R) —> L an embedding. If the image of ip is 
dense when projected to L/A, then for any bounded continuous rj 



^hm f r]{ip{^tUc.))diy{a) = f r]{M)dfi{M), 



;i) 



R'i-i 



L/A 



where v is any absolutely continuous probability measure on Ua and fi is the Haar probability 
measure on L/A. 

Remark 2. In effect, the Theorem says that the unstable manifold Ua is equidistributed in the 
larger homogeneous space L/A provided the density assumption is satisfied. 

We set A^ = e*-"^"^^* and apply the Theorem [3] with L = Ly., A = Ay, and 

(^:M^(l,«\...,«")-i(M,0,...,0)(l,i;\...,i;"). 

This ensures density after projecting to Ly j Ky by construction. We now use the following 
elementary Lemma to construct appropriate functions 77. 

Lemma 4. Let Nq = N U {0}. Then p2 : Nq x Nq — )■ Nq given by 



(x,y) ^ 



(Z/ + 1) 



is a bisection. 



By induction, there exists a polynomial bijection between Nq and Nq for each n\ call it 

Pn- Set 

/i„(M, a;\ . . . , a;") = p„(^,i,^i(M, a;^), . . . , ^.n,.n(M, x^)) 

and apply Shah's Theorem to functions 

^a(N, y\ . . . , ?/") = X{h„(M,a,i,...,3;")=p4Ai,...,A")}((l, *^\ • • • , V")(N, y\ . . . , y")) 



for all nonnegative integers AK These functions are not continuous, but are indicators of nice 
sets. Using a standard approximation argument we can apply the Theorem to them as well. 
For M of the form 



M 



/jY-i/(d-i) 



V 



]\[-l/(d-l) 



\ 



Nj 



(\ 



\ 



\(Xx ■ ■ ■ (Xd-\ 1/ 



we recover the numbers of points in the n segments. In fact, for all M 
r/A((l, «\ . . . , «")-'(M, 0, . . . , 0)(1, «\ . . . , «")) = X{,„(m,.i,...,.")=p4A)}(M, M«\ 

and 

^,,,(M,Mv)= ^ xi---XiX(r,r+a)(M(m + v)), 

and for the particular choice of M above. 



M('m. + v) 



md_i/Ar-V(«(-i) 
\{mxOLx H V nid-iad-i + ^Nj 



Mv' 



Thus, TjAi'piM)) = 1 if and only if there are exactly A^, . . . ,A"' visits to the segments around 
^^, . . . , ^" for the given ck and A^. Hence the form of the limiting distribution is 



pl'U,x/.(^^ • • • , ^") = h./aAUAm, x\ 



X 



= (Ai,...,A")}. 



(2) 

n 



3 Large d Limit 



In this section we consider the large d limit of the distributions from the previous sections 
and prove Theorem 2. Before proving the theorem, we will need basic information about the 
Poisson distribution. 

Poisson distribution with parameter a weighs each non-negative integer k with weight 
e~^a^ /k\ We will denote Poisson distribution with parameter a by Poisa. Its moments have 
the form 



oo 



n 



fe=0 ^ ^ fc=l 



S{n,k)a^ 



where S{n^ k) is the Stirling number of the second kind. As can be easily seen from the above 
equality, the Stirling number is the number of partitions of a set of n elements into k nonempty 
sets. The first few moments of the Poisson distribution are o", a"^ + a^ a^ -\- Sc^ + a. These 
correspond to partitions {!}; {12}, {11}; {123}, {112}, {121}, {211}, {111}. 



To further study the hmiting distributions we have obtained, we will need the following 
generalization of a proposition of Marklof from pQ which goes back to a theorem of Rogers [5] . 
Let Gr(n, /) = 0{n)/{0{l) x 0{n — /)) denote the Grassmannian of /-planes in R"; we assume 
that the /-planes are embedded in R" with respect to the standard basis. Let Gr(?T,,/)(Q) = 
{n G Gr(n, /) | n C R" is defined over Q} = {tt G Gr(ra, /) | vr fl Z" is a lattice in n}. For 
TT G Gr(r7,, /)(Q), we write covoIttz for the covolume of the lattice ttz = tt fl Z" in tt. We also 
set G = SL{d, R), r = SL{d, Z), and fix /i to be the Haar probability measure on G/T. 

Theorem 5. Let F: (R'^)®'' -^ H be a bounded piecewise continuous function with compact 
support. Let f : G/T — t- R fee defined by 

/(M) = J2 F{Mm\...,Mm'') 

with r < d a positive integer. Then, the first moment of f is given by the following expression: 

[ f{M)df,{M) = J2 E /^(^)r^r^' (3) 

J f^ crT1^^o^ ^ COVol TTz f 

where n' G Gr(r(i, /(i)(Q) is the image of tt under the embedding 

/"I 7"\/l 1 V V\ 

I nf nn 1 i L ( ly HT IT IT 1 

yjj , . . . , a. ) I r yd^ , . . . , j. , . . . , a. , . . . , a. ) 

and the measure dx is the Lebesgue measure on tt' that should be interpreted as the delta 
measure at the origin when 1 = 0. 

Remark 3. If in the sum over m^ we omit the terms where any of the m^ are 0, then in the sum 
over IT G Gr(r, /)(Q) we omit planes that are generated by subsets of the standard basis. This follows 
from the fact that such subsets of (Z )^ are SL((i, Z)-invariant. 

Lemma 6. With notation as in the Theorem, we have 

I Y^ F(Mm\...,Mm")d/i(M) = j F{x\ . . . ,x'')dxK . .dx\ 

G/r rn\...,rn';eZ'i ^,^-^d 

hnearly maep. 



Proof. First note that the integral is well-defined since linearly independent sets of vectors 
are preserved by F. Further renormalize /x so that fi{G/r) = 11^=2 C(^) fo^' d ^ 2 and write 
dfi{M)/fi{G/r) in the integral; this normalization will be useful later. Write 

I Xii . . . Xid\ 



M 



gG/F. 



\Xdl ■ ■ ■ XddJ 



Then for 1 < r < d we have 



M 



f a;ii . . 


Xir 


Orx(d-r) 




Xrl ■ ■ 


Xrr 






Xr+1,1 ■ ■ 


Xr+l,r 


dei"^{Xij)i^j^r .. 


. 


Xr+2,1 ■ ■ 


Xj,+2,r 


1 




\ Xdl 


Xdr 





1 



\ 



/ 





Zii . . 


■ Zi^d-r \ 


iCl^xr 








Zrl .. 


^r,d—r 




yii ■ 


yi,d-r 


0((i-r)xr 








yd-r,l ■ 


■ yd~r,d-r / 



where {i/ij) G SL{d — r, R). In these coordinates 



«i/i = J]^ dxij Y\ dzij ■ 6{l- det{yij)) J]^ dyij. 

i,j^d~r 



(4) 



i^d i^r 

j^d—r j^d—r 



The last factor is the Haar measure on SL{d — r, R) normahzed to ({2) . . . ({d — r) (or simply 
1 in case d — r = 1). 

For j = 1, . . . ,r let t-' = gcdm-'. Writing m^ /P for a column vector, we can find N G 
SL{d, Z) such that 



m 

71 



!rLi=N 



fan ai2 ■ 

022 • 

■ 

■ 



air\ 

a2: 







NA. 



V ■■■ / 

A is a matrix with integer entries uniquely determined by the following conditions: 

• ttij, 1 ^ i ^ j are relatively prime for any fixed j G {1, . . . , r} (in particular, an = 1); 

The first condition is due to relative primality of m^ /P , and the second comes from applying 
row operations. Given an = 1, a22, • • • , cirr, the number of possible matrices A of this form is 



where ipk is the number-theoretic function defined by 



Mf)=f'{^- 



1 

pk 



for A; ^ 1 and y^o is identically 1; y^i is Euler totient function. The function ipk{n) counts the 
number of /c-tuples (ni, . . . , rik) G {0, . . . , n — 1}^ such that gcd(n, rii, . . . , n^) = 1. 
Let us compute the stabilizer of a fixed matrix A: 



r^ = {7 e r I 7^ = A} 



idrxr 



-r)xr 



z. 



rx{d—r) 



SL{d-r,Z) 



Thus we get 



G/r "*' '• i- 



^^ ' ' ti,.,t'-=ig"^P Ner/FA 







MN 



/ai2\ 

022 



(air\ 



MN 



/i(G/r) 



E / ^(M 



voy 

/aii\ /ai2\ 

022 






G/Fa 







M 



)rf/i(M) 
/air\ 



,...,M 



VO/ 



V: / 



)rf/i(M). (5) 



Using dl]) to change the measure we get 



C(d)...C(d 



Id-r+i) y~! y~l / -^(^^33^ t^ai2aj^ + t^a22x'^-, • • • , fair^;^ H h fa„x'')dx^ . . . dx''. 



ti A 



(R'*)'- 



Now we do a linear change of variables and get 



(ti...r)'* Z^ a: 



C(rf)...C(rf-r'+i) Z^ (ti...r)'* Z^ a^2 ■ ■ ■ ""* 



a22) yr-l(arr) 



V=l a22v,arr = l 



■ f F{x\...,x'')dxK..dx\ (6) 



It is easy to see that 






^k{n) ^ C{d-k) 
n' ~ ad) ' 



whence the constant in front of the integral in ([S]) is 1, as desired. 
Proof of Theorem Rewrite the integral we are evaluating as 



D 



^^ ' ' J,-^ 1=0 rk(mi,...,m'-)=i 



Here we normalize /i like in the Lemma and the inner sum runs over those r-tuples of vectors 
whose R-span has dimension /. Since the set {rkm = /} is SL((i, Z)-invariant for each /, we 
can pass the sum over / through the integral sign. To prove the Theorem, it suffices to show 
that corresponding terms in the expression above and in (|3]) match for each /. Now observe 
that 

{rkm = 1} = M {r-tuples of vectors from (tt fl Z^Y with rank /}. 

7reGr(«,r)(Q) 

In fact, the sets whose union we are taking constitute an SL{d, Z)-invariant partition. Thus 
we need to parametrize linearly independent vectors in [ttHZ^Y for each tt. Let B: R' — )■ R^ 
be a linear map with image tt such that -B(Z') = tt n Z'". There can be many of these; any 
one will do. Using the standard basis, B = {hi) with 1 ^ i ^ / and 1 ^ j ^ r, and we obtain 
m^ = ^j 6^n*, where n* G Z*^ form a linearly independent set. Thus the integral becomes 



f,(G/r) 



E 

(//P 7reGr(«,r)(Q) 



J2 F{Mj2b]n\...,MY,bln^) 



Ti\...,n'gZ'* 
linearly indep. 



rf/i(M). 



The quantity in brackets is F-invariant, so the sum over vr can be interchanged with the 
integral. 

For each vr and B we can now apply the Lemma. It gives 



E / FiY,blx\...,Y,bW)dx'---dxK 



Since B(7}) = tt fl Z*", the Jacobian of B is the covolume of tt. The statement of the Theorem 
follows after a linear change of variables. 



D 



Proposition 7. For any S and r we have 



>(d) 

l,o-,H,r/H 



Poiso" 



as (i — 7- oo. 



Proof. From ([2]), all we need to show is that moments of fr^a are Poissonian for large d. First 
consider that case when ^ ^ Q. Without loss of generality we set r = 0. Then we need to find 



lim 

d—^oo 



MeG/r V I 



-1 k 



Yl iXi- ■■■■Xi-Xa)iMm + v) 

Tn6Z'*\{0} 



dfi{M)dv 



(7) 



for k = 0, 1, . . . Taking the integral over v inside the sum, we clear the way for Theorem [5] 
applied to 

F{xi, ...,Xd) = Gi{xi) . . . Gi{xd-i)Ga{xd) 



9 



where 

1 



Gt{z) = / xt{zi + y) ...xt{zk + y)dy. 
For any plane vr' as in the Theorem, we have that 



covoItt / f Gi{xd)dxd 

tt' 

It is elementary to see that the quantity raised to the power d is at most one: the numerator is 
the volume of vr "lying" inside the "crystal" shape, and covol ttz is the volume of a fundamental 
domain. To wit, consider first the case when Gt{z) is replaced by the indicator of [0, 1]'^. Since 
vertices of any fundamental domain for ttz have integer coordinates, it can completely cover 
the part of the plane inside the cube. Furthermore, the quantity in parentheses can equal one 
only when vr fl [0, l]'^ constitutes a fundamental domain for ttz- This means that there exists 
a Z-basis {ei}^ for vtz such that 

• Ci G {0,l}^ l^i^k; 

• Ci + ee e {0,l}^ lf^i,i'f^k. 

Hence two distinct e^, Cj' cannot take on the value 1 in the same coordinate. The same 
argument extends to other cubes of the form [— s, 1 — s]^ for s G [0, 1] and so, too, for the 
original Gt{z) as it is an average over cubes of this kind. 

The above argument shows that the limit as c/ — ;■ cxd exists for each moment and that 
rate of convergence is exponential. To understand this limit, we focus on the terms with 
J G{xi)dxi/ covolir = 1. Since in ([7]) we omit the terms in which any of m' = 0, the 
only terms that survive after taking the limit are the ones with tt generated by Cj for which 
"^i^i Cj = (1, . . . , 1) (no zero coordinates). For planes vr of fixed dimension / the number of 
possibilities is the number of partitions of a set of k elements into / non-empty subsets, which 
is exactly S{k, I). Finally observing that the last factor in ([8]) is u'^™'^ = a'', we find that the 
k-th moment tends to 

k 

J2s{k,iy, 
1=1 

which is the corresponding moment of the Poisson distribution with parameter a. 

In the case when ^ = p/q G Q we can modify the above proof. The integral over v becomes 
a finite sum, and we let Gt{z) = - Y^^Zq Xt(^i + ^/q) ■ ■ ■Xt{zd + r/q); the variable r appears in 
an equation similar to (|H]) and doesn't enter the definition of Gt{z). The statements from the 
continuous version are true for this function as well (since it is also an average over cubes), 
and the proof is complete. 

D 

Generalizing this proposition we can obtain the statement of Theorem 2. 
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Proof of Theorem 2. What we need to show is that 

n k^ 



: GZd\{o} 

has a hmit as d — ?■ oo for every choice of fc^, . . . , fc-^. If Y^^^^~ = (Y^, . . . ,Y'^) is dis- 
tributed according to PI^^.^^,^, then this expression is nothing more than the moment of 
order (/c^, . . . , fc-'). 

Now we make a simphfying observation: we can assume that y = 1 for all j without loss of 
generality since taking all possible n and S and computing E Y[ Y^ produces all the moments 
E'lK^"')'^^- That is, duplicating the random variable Y^ k^ times allows us to assume that k^ 
is 1. So we need to analyze 

/ / Yl Ylxi-----Xi-X{r^,a^+MMm^ + v^)dfi{M)dV, 

M y4d m^eZd\{0}i=l 



which by Theorem [5] is 



E' 



dx 



n ^1 • • • X(r^r^+<T^)(a;^ + v^)dV. (9) 



^ (covoIttz, 

7reGr(r,0(Q)^g=d ;;/ ^ ' 1=^ 

Since m^ are non-zero, we exclude the "coordinate planes" as in Remark [31 this is denoted 
by the prime in the formula above. 

We need to account for planes vr G Gr(r, /)(Q) that will contribute in the limit d — )■ oo. 
By the argument from the previous proposition we have that 



/n 
Ylxi X(T^,T^+aJ){x^ + v^)dx ^ (covoIttz)'^. (10) 



Since S*^ is normalized to have measure 1, it suffices to study the integrand for fixed V G S*^. 
If we can find V and tt for which strict inequality is true in (fTOj) , then by continuity we have 
strict inequality for the integral over V^ G S'^ and thus conlude that vr doesn't contribute in 
the limit. We will do this for V^ = first. A plane that will contribute in the limit c? — )■ oo 
must satisfy the property that vr fl [0, 1]^ is a fundamental domain for vr n Z^ as in the previous 
proposition. For each of these planes we can try to find another V that gives strict inequality 
in f fTOj) . If y G TT, we are translating the cube along the plane and thus getting the same 
cross-sectional area. So suppose V eE\7i; this corresponds to cutting the cube with a plane 
parallel to vr. It is easy to see that for such planes the section will always have smaller area 
than the one for V ^ tt. Thus it must be the case that S C tt. 

To summarize, a plane vr contributes to the limit only if vrfl [0, 1]^' is a fundamental domain 
for TT n Z*" and H C vr. This means that P „ ~ has a limit as c? — > oo because all moments 
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exist. If we write ^ = {$,^, . . . , ^^, ^^, . . . , ^^, . . . , ^"^ , . . . , ^"^ ) reordering as necessary, then vr 
must a product of admissible planes for (^^, . . . , ^^), (^^, . . . , ^^), . . . , (^" , . . . , ^" ). Hence 
the moment will split as the product of moments over distinct C,-^ ■ Using this observation 
and the previous proposition we see that in the case of distinct C,-' in Theorem 2 the limiting 
distribution is the product of independent Poisson distributions. If ^^ = ^^ but (r-', r-' + a-') fl 
(r-^ ,T^ + a^ ) = 0, then such ^^ and ^^ behave as if they were unequal since the factor 

from flTOj) vanishes. It is evident that if {t\ t^ + a^) fl (r-^ , r^ + a^ ) y^ for some j,j', then 
the limiting distribution cannot be a product of independent distributions. This concludes 
the proof of Theorem 2. D 
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